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Two  portions  of  a synthetic  hydroxyapatite 
(HA),  Ca5 OH (P04 )3  , fully  characterized  by  x-ray 
infrared,  petrographic,  and  chemical  analyses, 
were  heated  at  1,000°C  in  air  and  steam  atmos« 
pheres,  respectively.  Solubility  isotherms  for 
these  two  samples  in  the  system  Ca  (OH)  2 P04- • 

Hs0  were  determined  in  the  pH  range  5 to  7 by 
equilibrating  the  solids  with  dilute  H3P04  solu- 
tions. Both  samples  of  HA  dissolved  stoichio- 
metrically.  The  activity  products  (Ca++)° 

(0H~)  (POJ)3  —obtained  by  an  adjustment  of  the 
measurements  (Ca  and  P concentrations  and  pH 
of  the  saturated  solutions)  subject  to  the 
conditions  of  electroneutrality , constancy  of 
the  activity  product,  and  stoichiometric  dis- 
solution—-were  3.7S  ± 0.4e  x IQ-5 8 for  the  steam- 
heated  HA  and  2.5!  ± 0o037  x 10”BB  for  the  air- 
heated  HA.  Allowance  was  made  in  the  calculations 
for  the  presence  of  the  ion  pairs  [CaHP04j  and 
[CaH2P04]+.  The  higher  solubility  product  for 
the  air-heated  HA  is  ascribed  either  ^o  a change 
in  the  heat  of  formation  brought/by  partial  de- 
hydration or  to  a state  of  fine  subdivision 
resulting  from  a disproportionation  reaction. 

The  solubility  product  constants  were  used  to 
calculate  the  points  of  intersection  (i.e., 
singular  points)  of  the  two  HA  solubility  iso- 
therms with  the  isotherms  of  CaHP04*2H20  and 
CaHP04 ? it  was  found  that  the  pH's  of  the  singu- 
lar points  for  the  air-heated  HA  were  a full 
unit  higher  than  those  of  the  steam-heated 
preparation.  Conditions  are  described  for  the 
precipitation  of  HA  crystals  suitable  for  solu- 
bility measurements . 


This  investigation  was  supported  in  part  by  Research 
Grants  DE-01819  and  DE-02659  to  the  American  Dental 
Association  from  the  National  Institute  for  Dental 
Research. 
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1.  INTRODUCTION 


Hydroxyapatite  (HA) , Ca50H(P04)3/  has  wide  biological 
and  industrial  importance.  Its  physicochemical  properties, 
in  particular  solubility  and  detailed  crystallographic 
structure,  are  directly  related  to  the  mineralisation  pro- 
cesses in  biological  systems.  The  various  reports  in  the 
literature  on  solubility  of  HA  are  contradictory.  Claris, 
on  the  basis  of  precipitation  and  dissolution  experiments, 
concluded  that  HA  has  a definite  solubility  product. 
Thermodynamic  functions  for  HA  (from  which  a solubility 
product  constant  can  be  calculated)  have  also  been 
reported2^!-^.  Other  report^6^0^  however,  claim  that 
systems  saturated  with  respect  to  HA  do  not  conform  to  the 
principle  of  a solubility  product  constant  based  on  the 
stoichiometry  of  the  salt.  Among  the  models  that  have 
been  advanced  to  explain  the  apparent  anomalous  behavior 
of  HA  in  aqueous  solutions,  those  advocating  formation  of 
complexes  on  the  surface  of  the  solid6^"  have  received 
considerable  attention.  These  models  implicitly  assume 
that  equilibrium  between  the  solid  and  the  liquid  phases 
is  not  achieved  vzith  the  crystalline  portion  of  the  HA, 
but  instead  with  surface  layers  having  quite  different 
chemical  composition  from  that  of  the  bulk  of  the  solid. 
Furthermore,  the  composition  of  such  layers  has  been  as- 
sumed to  be  constant  over  a wide  range  of  calcium,  phos- 
phorus, and  hydrogen  concentrations  in  the  liquid  phase. 

At  least  part  of  the  problem  in  the  interpretation  of 
data  on  solubility  of  HA  is  related  to  the  difficulty  in 
preparing  pure  samples  of  this  compound.  So  far,  the  best 
crystals  have  been  obtained  by  hydrothermal  methods,  but 
the  yields  are  small,  the  equipment  is  relatively  expen- 
sive and  the  product  may  contain  extraneous  phases  such 
as  3-calcium  pyrophosphate  and/or  e-tricalcium  phosphate. 
Direct  precipitation  methods  usually  yield  very  finely 
divided  products,  often  without  the  stoichiometric  compo- 
sition of  HA.  A particularly  serious  difficulty  is  con- 
tamination by  precipitation  of  the  more  acid  calcium 
phosphates,  among  which  octacalcium  phosphate  (OCP) , 

Ca8H2 (P04)s • 5HS  0 and  its  decomposition  products  are  the 
most  difficult  to  eliminate. 
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The  purposes  of  this  work  were  (a)  to  devise  a simple 
method  for  preparation  of  HA  that  would  yield  sufficiently 
large  and  pure  crystals  so  that  an  unambiguous  character- 
ization would  be  possible,  (b)  to  learn  whether  this  com- 
pound behaves  as  a normal  chemical  compound  with  unique 
thermodynamic  properties  when  equilibrated  with  aqueous 
solutions,  and  (c)  to  establish  its  solubility  isotherms  and 
solubility  product  constant  at  25°Ca 


2.  EXPERIMENTAL  METHODS  AND  PROCEDURES 

In  the  method  of  preparation  outlined  here,  the  possi- 
bility of  OCP  precipitation  is  minimized  by  the  use  of 
relatively  high  temperature  and  pH's.  However,  pH  values 
higher  than  10  are  avoided  because  they  seem  to  favor  the 
formation  of  very  finely  divided  precipitates.  Other  con- 
ditions are  (a)  slow  addition  of  reagents  to  reduce  con- 
tinuous nucleation,  (b)  high  ionic  strength  and  the  presence 
of  acetate  ions  (a  Ca  complexing  agent)  to  increase  the 
solubility  of  HA,  and  thus  to  favor  formation  of  larger 
crystals,  (c)  use  of  a large  reaction  vessel  to  perform  the 
sacrificing  yields,  and  (d)  special  precautions  to  avoid 
contamination  with  carbonate  from  the  atmosphere  or  reagents 
(ammonia  was  used  to  control  the  pH  instead  of  NaOH  or 
KOH  because  the  alkali  hydroxides  are  difficult  to  obtain 
free  of  carbonate  contamination) . 


2.1.  Preparation  of  HA 

Freshly  boiled  distilled  water  and  chemically  pure 
salts  were  used  to  prepare  2.5  liter  of  7.25  x 10“2  M 
(NH4)2HPG4  solution  and  2.5  liter  of  13.33  x 10“2  M 
Ca (N03 ) 2 °4H2 0 solution;  their  pH's  were  adjusted  with 
ammonia  gas  to  8.5-9.  These  solutions  were  stored  in 
pyrex  bottles  vented  through  ascarite  tubes. 
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The  reaction  vessel  was  a 22-liter , 3-neck,  pyrex 
flask*  A 4-foot  condenser  was  attached  to  one  neck? 
transfer  of  reagents  and  materials  was  done  through  the 
other  two  necks*  The  top  of  the  condenser  was  connected 
by  a glass  tube  to  a 6-liter,  3-neck  flask  containing  3 
liters  of  10  N NaOH  solution?  a check  valve  permitting  out- 
ward flow  of  gases  was  connected  to  one  neck  and  a tube, 
extending  into  the  NaOH  solution  and  acting  as  a relief 
valve,  entered  through  the  third  neck.  The  reactant  solu- 
tions were  introduced  into  the  large  vessel  by  a peristal- 
tic pump  at  a rate  of  100  ml/hr.  Ammonia  gas  was  passed 
through  an  ascarite  tower  and  then  bubbled  through  water 
before  entering  the  reaction  vessel.  Samples  were  with- 
drawn through  a tube  that  extended  to  the  bottom  of  the 
reaction  vessel  for  petrographic  examination  of  the  solid 
and  pH  determination  of  the  liquid. 

Initially,  9 liters  of  N ammonium  acetate  solution 
were  boiled  in  the  reaction  vessel  for  about  12  hours  to 
remove  C02  ? then  ammonia  gas  was  passed  in  until  the  pH 
reached  a value  between  8.5  and  9 , and  the  pumping  of  the 
phosphate  and  calcium  solutions  was  started,  maintaining 
vigorous  boiling.  The  pH  was  maintained  between  8*5  and 
9.5  (measured  at  room  temperature)  by  an  adequate  flow 
of  NHg  until  the  total  volume  of  the  solutions  was  pumped 
into  the  reaction  vessel;  then  the  flow  of  NRj  was  dis- 
continued but  boiling  was  continued  for  about  three  hours. 
The  solid  was  allowed  to  settle,  the  supernatant  liquid 
was  withdrawn,  about  9 liters  of  boiled  distilled  water 
were  introduced  and  the  mixture  was  boiled  for  about  two 
hours.  This  washing  operation  was  repeated  as  'was  deemed 
necessary  to  remove  the  soluble  salts  (usually  five  times) * 
Further  washing  was  done  by  decantation,  using  boiled 
distilled  water  and  short  settling  periods  so  as  to  discard 
as  much  as  possible  the  finest  fraction  of  the  precipitate. 
Precautions  to  avoid  contamination  by  atmospheric  CC2 
were  taken  during  all  these  operations.  The  suspension 
was  then  filtered  and  the  solid  was  dried  at  150 °C  for  48 
hours.  As  explained  in  section  4,  a portion  of  this  sample 
was  heated  at  lyOOO^C  in  air  and  another  portion  was  heated 
at  the  same  temperature  in  steam  at  one  atmosphere* 
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2.2.  Equilibrations 

Equilibrations  were  made  in  glass-stoppered  bottles, 
each  of  which  contained  5g  of  HA  and  120  ml  of  a phos- 
phoric acid  solution,  by  slowly  rotating  end-over-end  in 
a water  bath  at  25  ± 0 . 1°C.  After  two  weeks,  the  suspen- 
sion was  transferred  to  a thermostated  saturator,  described 
in  an  earlier  publicatioiT^<  in  which  the  pH  of  the  clear 
filtered  solution  was  measured.  Samples  of  this  solution 
were  taken  for  calcium  and  phosphorus  analyses.  Standard- 
ization of  the  pH-meter  (with  claimed  relative  accuracy 
of  ± 0.0037  pH  units)  was  made  with  certified  NBS  buffer 
standards^'. 

Phosphoric  acid  solutions  used  in  the  equilibrations 
were  prepared  from  doubly  crystallized  2^P04*H20  and 
freshly  boiled  conductance  water.  The  molarities  of  the 
final  solutions  were  checked  by  chemical  analyses. 

2.3.  Calcium  and  Phosphorus  Determinations 

Calcium  was  determined  colorimetrically  by  the  in- 
direct technique  described  by  Banergee  _et  al^.  Phos- 
phorus was  determined  using  the  vanadomolybdate  reagent 
of  Brabson  et  a 3^3-.  in  both  analyses,  samples  were  com- 
pared directly  with  blanks  and  standards  in  a double-beam 
spectrophotometer . 

The  analytical  errors  estimated  for  the  measurements 
reported  here  were:  pH,  ± 0.009  units;  calcium  ± 2.5  percent? 

and  phosphorus,  ± 1.5  percent  of  the  amount  analyzed. 

3.  CHARACTERIZATION  OF  THE  SOLID 

The  HA  was  characterized  by  petrographic,  gas-adsorption, 
chemical,  infrared,  and  x-ray  methods. 


Under  the  petrographic  microscope,  the  crystals 
appeared  as  masses  of  birefringent , irregularly  shaped 
eguant  particles  up  to  5 ^ in  diameter.  Present,  also, 
was  a minor  fraction  of  thin  needles  up  to  10  m*  in  length. 
The  mean  index  of  refraction  (white  radiation,  25 cC)  was 
1.636  ± .002.  The  material  was  petrographicaXly  homo- 
geneous, both  the  equant  crystals  and  the  needles  are 
characteristic  of  HA.  The  mean  index  is  lower  than  the 
indexes  usually  given  for  HA  prepared  hydro thermally, 

Ne  = 1.644  and  N(JJ  = X.65137  , but  this  appears  to  be  a 
characteristic  of  precipitated  HA®  The  treatments  at  3/)00°C 
did  not  change  the  mean  indexes  significantly. 

The  surface  area  determined  by  the  BET  procedure*, 
using  nitrogen  as  the  adsorbate,  yielded  a surface  area 
of  6 ms/g.  Although  this  area  is  5 to  10  times  larger 
than  that  expected  from  microscope  examination,  the  agree- 
ment is  reasonable  considering  the  uncertainties  inherent 
to  the  two  methods. 

Chemical  analyses  yielded  18.22  ± .27  percent  P and 
39.63  ± .99  percent  Ca  as  compared  to  the  theoretical 
values  18.53  percent  P and  39.84  percent  Ca  for  HA. 


*Measurement  made  by  Dr.  W.  V.  Loebenstein,  MBS 
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The  infrared  spectrum  displayed  all  the  bands  charac- 
teristic of  well-crystallized  HA^.  Small  bands  in  the 
regions  1410  and  1460  cm"1  indicated  slight  contamination 
with  carbonate  and  moisture,  respectively.  Analysis  by 
the  Warburgh  method  yielded  0.04  percent  COs  . 

The  x-ray  diffraction  pattern  was  sharp  and  character- 
istic of  that  for  pure  No  extraneous  peaks  were 

found. 


Although  the  above  methods  to  characterize  the  HA  did 
not  reveal  the  presence  of  extraneous  phases,  when  a 
sample  of  the  HA  was  heated  to  1,000°C  for  24  hours,  the 
infrared  spectrum  showed  small  peaks  in  the  region  920  to 
1,200  cm-*1.  These  peaks  indicate  the  presence  of 
(3-Cag  (P04)2  14,  which  presumably  formed  on  heating,  from  a 
contaminant,  an  acid  calcium  phosphate. 

A portion  of  the  HA,  therefore,  was  heated  at  1,000°C 
for  six  days  in  room  atmosphere  to  anneal  the  crystals, 
to  drive  off  carbon  dioxide  and  moisture,  and  to  convert 
the  less  basic  calcium  phosphate  into  $-Ca3  (P04)2.  Following 
the  heat  treatment,  the  sample  was  extracted  with  2 x 10"3  M 
Hg  P04  to  remove  the  B-Cag  (P04)2.  Approximately  10  percent 
of  the  sample  was  dissolved  by  this  treatment. 

Another  portion  of  the  HA  was  subjected  to  the  same 
heat  treatment  but  in  the  presence  of  one  atmosphere  of 
water  vapor,  and  then  extracted  with  the  P04  solution. 

The  infrared  spectra  of  the  heated  samples  were  very 
similar,  both  corresponding  to  pure  HA^>  but  the  height  of 
the  peak  at  3,570  cm"1  was  considerably  greater  for  the 
steam-heated  sample,  indicating  either  a more  uniform  en- 
vironment for  the  0H“  group  in  this  sample  or  actually  a 
higher  content  of  OH”  in  the  crystals. 
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Another  sample  of  HA,  referred  to  here  as  "as  prepared1' 
was  synthesized  according  to  the  method  described  above, 
but  it  was  not  subjected  to  high  temperature  treatment. 


4.  CALCULATION  OF  SOLUBILITY  PRODUCT  CONSTANTS 
The  solubility  product  constants. 


where  parentheses  denote  ionic  activities,  were  calculated 
by  a procedure  very  similar  to  the  one  described  elsewher£?< 
The  values  are  listed  in  tables  1,  2 and  5,  However,  it 
was  recognized  that  the  best  values  to  be  derived  from  such 
data  must  take  into  account  more  completely  the  propagation 
of  errors  and  the  proper  weight  of  the  measurements.  A 
more  refined  method  was  devised,  therefore,  for  these  cal- 
culations which  entails  the  simultaneous  adjustment  of 
the  weighted  observables— -i.e. , calcium  and  phosphorus 
concentrations  and  pH's  in  the  various  equilibrated  solu- 
tions plus  the  initial  concentration  of  Hg  P04— and  the 
determination  of  a parameter  which  yields  best  estimates 
of  Ksp.  The  main  features  of  the  method  follow  the  gen- 
eral least-squares  adjustment  discussed  by  Deming3^  and 
they  are  similar  to  those  recently  presented  by  Wentworth 
et  al3^  for  the  calculation  of  complexing  constants  using 
spectrophotometric  data.  However,  for  the  present  work, 
two  or  three  condition  functions  were  used  instead  of  one. 

A detailed  account  of  the  method  will  be  given  elsewhere3^. 
It  suffices  here  to  give  the  three  condition  functions  used 
in  this  investigation. 


Ksp  = (Ca++ ) 5 (OH-) (P04*)3 
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in  which  M and  P are  the  total  molar  concentrations  of  Ca 
and  P in  solution,  respectively 7 (H+)  is  the  activity  of 

the  hydrogen  ion?  f is  the  molar  activity  coefficient  of 
the  subscripted  species?  Kw  and  Kg  are  the  dissociation 
constant  of  water  and  the  third  ionization  constant  of 
Hg P04 , respectively 7 P0  is  the  initial  molar  concentration 
of  Hg  PO4  7 S is  the  parameter  to  be  calculated,  its  value 
being  the  best  estimate  of  Ksp?  the  fraction  5/3  is  the 
theoretical  Ca/P  ratio  in  HA?  T is  the  total  molar  concen- 
tration of  Ca  or  P present  as  the  ion  pairs  [CaHP04]°  and 
[CaHsPQ4]+  V and  it  is  defined,  as  are  the  quantities  R, 

Q1  and  Qg , by  the  expressions, 

T = (1/2  M + 1/2  P + L)  -7(1/2  M + 1/2  P + L)a  - (MP)  [4] 
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In  equations  [4-7],  K:  and  K2  are  the  first  and  second 

ionization  constants  of  Eq  P04 , respectively,  and  Kx  and 
Ky  are  the  association  constants  for  the  ion-pairs  [CaHP04]° 
and  [CaH2P04  ] + . A list  of  the  numerical  values  for  the 
various  constants  used  in  this  investigation  is  given  in 
the  Appendix. 
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The  three  condition  functions,  equations  [1]  to  [3], 
are  based  on  electroneutrality,  the  solubility  product  for 
HA,  and  congruency  in  the  dissolution  of  the  solid,  respec- 
tively. These  three  conditions  were  used  to  obtain  the 
best  estimates  of  Ksp  for  the  preheated  hydroxyapatites  (in 
air  and  steam) . Condition  [3] , defining  stoichiometric 
dissolution,  was  excluded  for  the  calculations  with  HA  "as 
prepared"  because  phosphate  was  adsorbed  in  some  of  these 
equilibrations,  as  explained  in  section  5.1. 

Equation  [4]  was  obtained  from  the  simultaneous  solu- 
tion of  equations  [8]  and  [9]  which  describe  the  equilibria 
involved  in  the  formation  of  ion  pairs.  Calling  X the 
concentration  of  [CaHP04]°  and  Y that  of  [CaH3F04]+, 


X = 


Kx  [M-t]  [P— t 3 f c a + 
02  f x 


[8] 


Y = 


Ky  [M-t]  [P-T]fc++(H+) 
02KSfy  " 


[9] 


in  which  fx  and  fy  are  the  activity  coefficients  for  the 
species  [CaHP04]°  and  [CaH2P04  ] + , respectively. 


by 


The  ionic  strength,  p,,  was  defined  in  these  calculations 


\i  = 2 [M-t]  + 1/2 


in  which 


(H+ ) 


+ 


K. 


(K+)f0  H- 


+ LT 


(P— t)  Qq 
Qo 


[10] 


03  = 


(H+ ) 


K2  fH2  P04 


+ 


+ 


£h  pd“  (H  )fp'o^ 


[11] 
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The  values  for  the  ionic  activity  coefficients  were 
calculated  by  the  use  of  the  Debye-Hiickel  theory  using 
the  distances  of  the  closest  approach  given  in  the  Appendix. 
The  iterative  procedure,  starting  with  unitary  activity 
coefficients,  reported  in  previous  publication was 
used  for  these  calculations.  The  activity  coefficient  fx 
was  taken  to  be  unity,  and  it  was  assumed  that  fy  = 

fH2>i04“* 


5.  RESULTS 

5.1.  Solubility  of  HA 

The  experimental  results  for  the  steam-heated  and 
air-heated  preparations  are  given  in  Tables  1 and  2, 
respectively,  along  with  the  values  of  Ksp  calculated 
from  these  results  using  the  activity  product  that  appears 
as  the  leading  term  in  equation  [2].  The  spread  in  Ksp 
in  each  table  is  acceptably  low  even  though  the  average 
Ksp  for  the  steam-heated  material.  Table  1,  is  three  orders 
of  magnitude  less  than  that  for  the  air-heated  material. 
Table  2.  The  other  tests  that  can  be  applied  to  these 
data— apparent  electroneutrality  unbalance,  stoichiometry 
of  the  dissolution  reaction,  and  smoothness  of  the  plots 
of  Ca  or  P vs.  pH— also  show  deviations  which  are  within 
the  anticipated  experimental  errors,  thus  vouching  for 
the  validity  of  the  results.  In  Tables  3 and  4 are  given 
the  adjusted  values  for  the  three  observables,  Ca,  P and 
pH,  and  the  adjusted  values  of  Ksp,  electroneutrality 
and  the  ratio  Ca/P  of  the  dissolution  reactions.  It  is 
clear  that  the  adjustments  in  the  Ca,  P,  and  pH  values  are 
within  their  anticipated  experimental  errors.  The  fact 
that  these  adjustments  are  small  is  further  confirmation 
of  the  internal  consistency  of  the  data.  We  believe  that 
the  best  estimate  of  the  two  Ksp 1 s obtained  by  the  adjust- 
ment procedure  (Tables  3 and  4) , which  weights  the  experi- 
mental measurements  according  to  their  errors,  constitute 
better  values  for  the  solubility  products  than  the  simple 
unweighted  mean  values,  given  in  Tables  1 and  2. 
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The  large  difference  between  the  two  Ksp  values 
(Tables  3 and  4)  probably  relates  to  structural  changes 
induced  by  the  heat  treatments,  as  discussed  later. 

Also  listed  in  Tables  3 and  4 are  the  values  for 
ionic  strength  and  the  fraction,  r,  of  Ca  present  in  the 
form  of  ion  pairs  given  by  t/M.  The  gradual  increase  in 
r with  decreasing  pH  relates  to  the  increase  in  Ca  and  P 
concentrations,  and  are  mostly  the  result  of  an  increase 
in  the  concentration  of  the  ion-pair  [CaH2P04]+. 

Results  of  the  solubility  experiments  for  the  11  as 
prepared"  HA  are  shown  in  Table  5.  The  large  variation 
in  the  Ksp  values  in  Table  5 is  believed  to  relate  to  the 
presence  of  other  calcium  phosphates  less  basic  than  HA? 
evidence  for  this  type  of  contamination  is  presented  in 
Section  4.  The  data  also  suggest  the  presence  of  an  anion 
foreign  to  the  system  Ca (OH) 2-H3P04-Hs0  as  evidenced  by  an 
apparent  positive  electroneutrality  unbalance  which  is  far 
beyond  experimental  error  for  the  three  treatments  with 
pH's  above  6.  This  may  have  been  caused  by_ carbonate 
since  the  infrared  spectrum  showed  small  COg  peaks. 

Interpretation  of  the  data  in  Table  5 is  further 
complicated  by  sorption  phenomena?  for  example,  the 
phosphorus  concentrations  in  the  solutions  for  the  last 
three  entries  of  Table  5 are  less  than  the  initial  H3P04 
concentrations.  A systematic  study  of  this  type  of  surface 
phenomenon  with  finely  divided  HA  has  been  reported  else- 
where^£~.  Although  some  of  the  Ksp  values  in  Table  5 
are  in  the  range  of  those  given  in  Table  1 for  the  steam- 
heated  HA,  the  three  at  the  higher  pH's  are  abnormally 
large.  Furthermore,  the  calculated  adjustments  of  the 
measurements,  not  shown  here,  were  so  large  as  to  invali- 
date the  least  squares  procedure.  As  a result,  little 
weight  is  given  to  the  set  of  data  in  Table  5.  They  give 
some  confirmation  in  the  low  pH  range  of  the  results 
obtained  for  the  steam-heated  HA  and  they  illustrate  how 
the  solutions  can  be  made  supersaturated  in  the  high  pH 
range  by  the  presence  of  calcium  phosphates  less  basic 
than  HA  and/or  contamination  with  carbonates. 
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The  phosphorus  concentrations  of  the  solutions 
saturated  with  respect  to  the  preheated  HA's,  Tables  1 
and  2,  are  plotted  against  their  corresponding  pH's  in 
Figure  1.  Similar  plots  for  solutions  saturated  with 
respect  to  CaHP04*2H.,0  (DCPDy^^-iiS-  and  CaHP04  (DCPA]'*^" 
at  25 °C  are  also  shown  in  the  same  figure*  The  HA 
curves  were  extrapolated  (dotted  lines)  to  meet  the 
intersection  points  A,  B#  c,  and  D in  Figure  1.  These 
four  points  correspond  to  singular  points  between  the 
preheated  HA  9s  and  DGPD  or  DCPA.  The  chemical  compo- 
sitions of  the  solutions  at  the  singular  points  were 
calculated  by  solving  equation  [12]  for  (H+)  through  a 
method  of  successive  approximations , 


2 (H*)  sKnf,K  (H+ ) RK.K,,.  _ K* = 

f«++  f fH+  + 1-R  " tiff*  (H+)fQH- 

in  which  Koc  represents  the  solubility  product  constant 
for  either  DCPD  or  DCPA,  as  the  case  may  be,  and  R and 
K are  defined  by 


R 


1 


1 + KxKgfH2PO-4 

Ky (H+ ) 


[13] 


~ Ksp  1 I 
.KfKJKlc J 


[14] 


The  value  for  (H+)  obtained  from  equation  [12]  was  then 
used  to  obtain  the  numerical  value  of  t since,  for  satura- 
tion with  respect  to  either  DCPD  or  DCPA, 


K,Koc 

" 1-H 


[15] 
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The  concentrations  of  calcium  and  phosphorus  are  obtained 
from  equations  [16]  and  [17],  respectively 0 


M = 


(H+) 


KDcK 


-p  + + 
-LCa 


+ 


[16] 


P 


q5 

(H+ ) 2 K 


[17] 


The  chemical  compositions  of  the  solutions  at  the  four 
singular  points,  calculated  in  the  foregoing  fashion,  are 
given  in  Table  6.  Considerably  higher  concentrations  of 
Ca  and  P (and  lower  pH's)  are  found  at  the  singular  points 
of  the  preheated  HA ' s with  DCPD  than  at  the  corresponding 
singular  points  with  DCPA;  the  same  pattern  is  observed 
for  the  fractions  of  Ca  present  as  ion  pairs  (r  in  the 
table) . These  results  are  a consequence  of  the  higher 
solubility  of  DCPD1^^-2-  relative  to  the  anhydrous  salt  at 
25°C The  most  striking  differences  in  the  compositions 
at  the  singular  points,  however,  are  brought  about  by  the 
solubility  of  the  HA  itself.  Thus,  a full  pH  unit  dif- 
ference exists  between  the  singular  points  of  air-heated 
and  steam-heated  HA  with  either  of  the  two  dicalcium  phos- 
phates; changes  in  calcium  and  phosphorus  concentrations 
by  factors  of  about  3.5  and  4,  respectively,  parallel 
this  pH  difference.  If,  as  it  is  apparent  from  the  present 
results,  the  solubility  of  HA  depends  upon  the  previous 
history  of  the  sample,  it  is  nearly  impossible  to  predict 
the  composition  of  singular  points  between  a precipitated 
HA  and  other  calcium  phosphates  without  fully  character- 
izing the  solubility  properties  of  the  HA  involved. 
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6.  DISaJSSION 
6.1.  Solubility  of  HA 


The  solubility  data  given  in  Tables  1 and  2 for  the 
two  preheated  HA  samples  are  consistent  with  stoichiometric 
dissolution  represented  by  the  reaction  Ca50H(P04)3 ^5Ca++ 

+ 0H'~  + 3P 04  = . It  is  believed  that  saturation  conditions 
were  obtained  in  the  two-week  equilibrations  as  evidenced 
by  the  smooth  curves  obtained  in  the  plots  of  calcium 
concentration  versus  pH  (Figure  1} , calcium  concentration 
versus  phosphorus  concentration  (not  shown  in  this  paper) , 
and  a reasonable  constancy  in  the  ionic  activity  product 
for  HA.  The  success  of  the  adjustment  procedure  described 
in  section  5 is  additional  evidence  that  the  parameters  S, 
Tables  3 and  4,  are  valid  solubility  product  constants 
for  the  preheated  hydroxyapatites. 

The  higher  solubility  product  constant  obtained  for 
the  HA  heated  in  air,  relative  to  the  product  obtained 
for  the  sample  heated  in  steam,  is  most  probably  due  to 
structural  changes  induced  by  the  thermal  treatment.  Such 
changes  might  relate  in  part  to  the  statistically  disordered 
structure  of  HA  described  by  Kay,  et  al^  originating  from 
small  displacements  of  the  OH-  groups  from  the  mirror 
planes  in  the  space  group  P63/m«  If  this  were  purely  an 
entropy  effect,  such  disorder  would  tend  to  stabilize  the 
lattice;  the  steam-heated  sample,  therefore,  would  have 
to  be  the  disordered  one.  However,  it  is  difficult  to 
see  why  heating  HA  in  air  would  induce  more  order  in  the 
OH"  groups  than  would  be  obtained  by  heating  in  steam. 

A more  likely  explanation  is  that  heating  in  air  results 
in  loss  of  OH"  groups  by  partial  dehydration  of  the  HA. 
Indeed,  the  infrared  spectra  tend  to  substantiate  this 
alternative.  Whereas  the  half  width  of  the  peak  at 
3570  cm”  corresponding  to  the  OH"  stretching  mode,  was 
essentially  the  same  for  both  preheated  hydroxyapatites, 
its  height  for  the  steam-heated  HA  was  almost  twice  that 
of  the  air-heated  HA. 
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Although  the  weight  loss  in  our  sample  was  not 
measured,  slight  weight  losses  have  been  reported^  in 
the  range  3/300 °C-3/100 °C  for  synthetic  HA  samples  heated 
gradually  from  25°  to  3/200  °C;  these  losses  occurred 
without  alteration  in  their  x-ray  patterns 0 Also,  it 
has  been  reported^-3"  that  at  3/200  °C  in  air,  dispropor- 
tionation of  HA  into  tricalcium  phosphate  and  tetracalcium 
phosphate  occurs  according  to  the  reaction, 

2Ca5  OH  (P04  )3  - 2 cc-Ca,  (P04)2  + Ca40(P04)2  + Hs  0 [18] 

The  x-ray  powder  diffraction  pattern  of  the  air-heated 
sample  was  the  same  as  that  of  HA,  but  it  is  impossible 
to  rule  out  a partial  disproportionation. 

In  view  of  the  foregoing  considerations,  there  are 
two  explanations  for  the  higher  solubility  of  the  air- 
heated  HA:  (1)  Reversal  of  the  disproportionation 

reaction  [18]  when  the  sample  is  returned  to  an  aqueous 
system  might  yield  very  finely  divided  HA  with  an  unusually 
high  solubility  because  of  surface  energy  contribution, 
or  (2)  loss  of  OH"  ions  by  dehydration  during  heating  of 
the  solid,  without  destruction  of  the  HA  lattice,  might 
yield  crystals  with  higher  free  energy.  Evidently, 
the  process  responsible  for  the  increased  solubility  is 
not  readily  reversible,  since  otherwise  the  values  for  the 
Ksp  of  the  two  preparations  w7ould  have  been  the  same.  There- 
fore, the  equilibria  attained  with  the  air-heated  HA  must 
have  been  metastable  relative  to  the  more  stable  form 
of  HA.  This  would  require  that  the  rate  of  dissolution  of 
the  air-heated  HA  was  higher  than  the  rate  of  precipitation 
of  the  more  stable  salt. 
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The  Ksp  values  for  the  two  kinds  of  preheated  apatites 
correspond  to  a difference  of  about  4 Kcal  mol-1  in  the 
standard  free  energies  of  formation.  If  it  were  assumed 
that  the  standard  heat  of  formation  is  the  same  for  both 
compounds,  the  observed  difference  in  free  energy  would 
correspond  to  about  14  cal  deg”1  mol-1  less  entropy  in 
the  air -heated  apatite  than  in  the  steam-heated  salt. 

This  difference  appears  too  large  to  be  associated  with 
a change  in  entropy  resulting  from  vacancies  or  ordering. 
For  this  reason  it  is  believed  that  the  difference  in 
solubility  is  mostly  related  to  a difference  in  their 
heat  contents. 


6.2.  Singular  Points 

A system  of  two  solid  phases  in  equilibrium  with 
an  aqueous  phase  in  the  ternary  system  Ca (OH) 2 -Hg  P04— H2  0 
is  characterized  by  zeto  degrees  of  freedom  (under 
conditions  of  constant  temperature  and  pressure) . Such 
singular  points  are  important  in  solubility  diagrams 
because  they  correspond  to  points  where  there  is  a change 
in  the  relative  stability  of  the  two  solid  phases.  Thus, 
with  reference  to  Table  6,  in  saturated  solutions  with 
pH's  below  4C33,  DCPD  is  more  stable  than  the  steam-heated 
HA;  in  solutions  with  pH's  above  4.33,  steam-heated  HA 
is  more  stable  than  DCPD.  The  greater  solubility  of  the 
air-heated  HA  has  the  effect  of  increasing  the  pH  of  the 
singular  points  with  the  dicalcium  phosphates  by  about 
one  pH  unit  relative  to  those  with  steam-heated  HA.  If 
an  apatite  comparable  to  the  air-heated  HA  were  formed 
in  oral  environments,  its  greater  solubility  would  extend 
the  stability  ranges  of  the  two  dicalcium  phosphates  to 
conditions  which  occur  in  saliva  in  vivo. 

In  the  concentration  ranges  covered  by  the  present 
experiments,  the  degree  of  ion-pair  formation  varied  from 
about  1„4%  (at  the  highest  pH)  to  about  3.4%  (at  the  lowest 
pH)  of  the  total  calcium  concentration 0 In  view  of  the 
uncertainty,  both  theoretical  and  experimental  in  calcu- 
lations for  non-symetr ical  electrolytes,  the  formation  of 
ion-pairs  in  the  experimental  segments  of  the  solubility 
isotherms  for  HA  is  not  considered  significant.  However, 
in  more  acid  solutions,  such  as  those  of  the  calculated 
singular  points,  the  ion-pair  formation  does  become 
significant. 
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6.3.  Stoichiometry  and  Potential  Diagrams 

The  stoichiometry  of  a solid  calcium  phosphate 
equilibrated  with  aqueous  solutions  is  easily  revealed 
by  the  use  of  "potential  diagrams These  are 
solubility  diagrams  in  which  the  coordinates  are  functions 
of  the  chemical  potentials  of  the  components , Ca(0H)2  and 
H,P04,  in  solution0  For  example,  if  the  quantity  - log 
(Ca^  + ) (OH-) 2 is  plotted  against  - log  (H+)3(P0j),  the 
saturation  condition  is  defined  by  a straight  line  and 
the  absolute  value  of  its  slope  corresponds  to  the  reci- 
procal of  the  Ca/P  ratio  in  the  solid.  This  kind  of 
diagram  is  illustrated  in  Figure  2 with  calculations 
based  on  the  data  for  the  two  preheated  apatites  (Tables 
1 and  2) . Similar  calculations  were  made  with  data  from 
other  investigators  and  these  are  plotted  in  the  same 
diagram.  The  slopes  of  the  regression  lines  for  the  air- 
heated  and  steam-heated  apatites  shown  in  Figure  2 corres- 
pond to  apparent  Ca/P  ratios  of  lc69  ± 0.02  and  1.68  ± 0o07, 
respectively o These  ratios  agree  with  the  theoretical 
ratio  for  HA,  1.67,  and  the  ratio,  1.68  ± 0*05  obtained 
by  chemical  analysis  of  the  untreated  solidc  In  view  of 
these  facts,  the  relatively  low  value  obtained  for  the 
Ca/P  ratio  in  the  dissolution  reaction,  1.57  ± 0o02, 
of  the  steam-heated  preparation.  Table  1,  must  be  ascribed 
to  experimental  errors.  It  is  to  be  noted  that  the  Ca/P 
ratio  calculated  from  dissolution  involves  a difference 
between  two  experimental  measurements. 


6.4.  Comparison  with  Reported  Data 

Clark^'  reported  a value  of  57.75  as  his  best  estimate 
of  pKsp  (the  negative  logarithm  of  Ksp)  of  HA  at  25 °C. 

(The  actual  value  he  reported,  115.5,  is  based  on  the  form- 
ula Ca7 o (OH) 2 (P04 ) s rather  than  the  formula  Ca50H(P04)3 
used  here.)  His  systems  consisted  of  HA  precipitated 
by  mixing  dilute  solutions  of  Ca(0H)2  and  HgP04  at  90°C; 
he  kept  the  resulting  suspensions  at  that  temperature 
for  120  hours  and  re-equilibrated  them  at  25 °C  for  90 
hours.  The  Ksp  corresponding  to  his  pKsp  value  of  57.75 
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is  1.78  x 10“58,  which  is  in  good  agreement  with  the 
value  3.73  x 10“58  reported  here  for  the  steam-heated 
HA.  The  individual  values  reported  by  Clark  range  from 
0.53  x 10  °8  to  7.59  x 10“58.  Even  though  Clark's  data, 
as  plotted  in  the  chemical  potential  diagram.  Figure  2, 
have  some  scatter,  it  is  apparent  that  they  approximate 
closely  the  same  line  as  our  data  for  the  steam-heated 
HA.  The  same  author  observed  that  25 °C  equilibrations 
made  with  HA  precipitated  at  40°C  resulted  in  large 
apparent  supersaturations  with  respect  to  the  values  for 
Ksp  given  above.  In  light  of  present  knowledge,  one 
would  expect  that  those  precipitates  contained  some  of 
the  more  acid  calcium  phosphate  which  are  more  soluble 
than  HA.  Thus,  the  more  reliable  of  Clark's  results 
agree  with  those  reported  here  for  steam-heated  HA.  We 
consider  it  significant  that  such  good  agreement  was 
obtained  between  materials  prepared  by  precipitation 
from  aqueous  systems  and  a material  subjected  to  anneal- 
ing at  1,000°C  in  a steam  atmosphere.  This  suggests  that 
the  Ksp  value  obtained  for  the  steam-heated  material  is 
characteristic  of  any  well-crystallized  HA. 

A solubility  product  constant  for  HA  was  calculated 
from  the  following  thermodynamic  quantities:  aH|  = - 

3,212  Kcal  mol”1  reported  by  Smirnova  et  al^-1-  fGr  HA 
which  had  also  been  heated  at  1,000°C  in  the  presence 
of  steam?  S°  = 186.6  cal  mol”1  deg""1  reported  by  Egan 
and  Wakefield4  for  another  HA  treated  in  a similar  way, 
and  standard  entropies  of  formation  of  the  elements  and 
standard  free  energies  of  formation  of  ions  reported 
in  NBS  Publications’2'---!^'.  The  product  thus  obtained  is 
of  the  order  !0~57. 

This  result,  too,  is  considered  to  be  in  good 
agreement  with  the  order  of  magnitude  given  here,  10*“58, 
for  the  steam-heated  apatite,  considering  the  uncertain- 
ties in  some  of  the  thermodynamic  quantities  used  and 
the  great  effect  that  small  variations  in  these  quanti- 
ties have  on  the  Ksp  calculated  in  this  fashion. 
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The  foregoing  results  are  not  in  agreement  with 
those  of  Levinskas^-,  Levinskas  and  Neuman^  and  Neuman 
and  Neuman^  who  concluded  that  synthetic  preparations 
of  HA  equilibrated  in  aqueous  solutions  do  not  obey 
the  solubility-product  constant  principle 0 One  of  the 
hydroxyapatites  (L-apatite)  most  extensively  studied 
by  these  authors  was  a commercial  preparation  of  fine 
particle  size,  68  m2/g.  The  large  surface  areas  of  their 
materials  could  introduce  difficulties  in  the  interpre- 
tation of  solubility  studies.  However,  neither  theoret- 
ically^-0 nor  experimentally^"'  do  sorption  phenomena 
invalidate  the  application  of  the  solubility  product  con- 
stant to  systems  equilibrated  with  HAC  It  should  be 
pointed  out  that  the  sample  of  HA  (L-apatite)  used  by 
Neuman  and  coworkers  had  a refractive  index  of  1.605s 
which  is  much  too  low  for  a well-crystallized  HA;  in 
addition,  their  preparation  was  contaminated  with  cal- 
cite,  contained  0o28  percent  sulfate  and  0.8  percent 
C0S  5 . No  history  is  given  for  the  synthesis  of  L-apatite, 
but  the  analytical  results  clearly  show  the  heterogeneity 
of  the  sample0  Under  these  circumstances,  reactions 
other  than  simple  dissolution  and  precipitation  of  HA  are 
to  be  expected,  and  their  involved  kinetics  is  a burden 
for  sound  interpretations  of  solubility  studies.  The 
adoption  of  the  ion-activity  product  (Ca+  + ) (KPO^)  to 
characterize  HA  systems,  as  suggested  by  Neuman  and  co*~ 
worker^?9  cannot  be  justified  on  theoretical  base  £3^ 

Our  recalculation  of  their  data^l-  plotted  in  Figure  2 
shows  that  the  points  are  widely  scattered,  do  not  define 
a single  line,  and  cannot  be  used  to  derive  a reliable 
Ksp  for  a solid  of  any  Ca/P  ratio. 
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Rootare,  Deitz,  and  Carpenter^  reported  results  on 
solubility  studies  done  with  several  commercial  and 
laboratory-prepared  samples . They  observed  wide  variations 
in  the  calculated  ion-activity  products  for  HA;  these  were 
explained  by  postulating  formation  of  a surface  complex 
according  to  the  reactions, 

Ca10 (P04)s  (OH) s + 6HsO^[Ca2 (HP04) (0H)2 ] + 2Ca++  + 2HPOJ  [19] 
Ca2 (HP04 ) (0H)-,^2Ca+  + + HPOj  + 20H”  [20] 

It  is  true,  as  they  claim,  that  the  ionic  product  calcu- 
lated according  to  equation  [20]  showed  less  variability 
(pK  values  from  24.6  to  2608)  than  those  calculated  for 
HA  on  the  basis  of  the  formula  unit  Cax 0 (OH) g (P04 ) s (pK 
from  108  to  120) 0 However,  the  mechanism  postulated  by 
Rootare  et  ad.  is  open  to  the  following  criticisms: 
a)  Although  the  composition  of  the  HA  surface  does  not 
have  to  conform  necessarily  to  the  stoichiometry  of  the 
bulk  solid,  there  is  no  a priori  reason  to  postulate  that 
certain  surface  compositions  will  prevent  equilibrium  with 
the  crystalline  HA  latticed;  b)  It  is  highly  improbable 
that  the  surface  of  HA  equilibrated  in  aqueous  solutions 
has  a definite  composition  independent  of  the  composition 
of  the  aqueous  mediuftr^";  c)  If  a definite  composition 
is  ascribed  to  the  "surface  complex"  and  the  solubility 
product  principle  is  applied  according  to  its  stoichiometry, 
it  is  mandatory  that  such  species  be  considered  as  another 
solid  phase  (i.e„,  has  constant  molar  free  energy  and 
stoichiometry^)  ; then,  in  the  light  of  equations  [19] 
and  [20] , the  system  would  have  zero  degrees  of  freedom 
(in  the  system  Ca (OH) 2 -H3P04-H3 O)  under  fixed  temperature 
and  pressure^  so  that  the  solution  composition  could 
not  vary  as  observed  by  them;  d)  It  can  be  easily  shown 
that  the  equilibrium  conditions  [19]  and  [20],  as  given 
by  Rootare  _et  al,  also  imply  equilibrium  with  respect 
to  the  HA  solid  phase;  e)  Of  the  12  samples  used  by 
Rootare  et  a].,  only  four  were  petrographically 
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homogeneous5 ; the  rest  were  contaminated  with  calcite 
varying  from  traces  to  about  one  percent,  with  sulfate 
varying  from  0.006  to  0.28  percent,  and  with  the  exception 
of  two,  all  the  samples  contained  carbonate  in  significant 
amounts;  only  the  four  samples  prepared  in  the  laboratory 
showed  a refractive  index  approaching  the  one  expected 
for  well-crystallized  HA;  one  of  the  latter,  P-20,  probably 
contained  less  basic  calcium  phosphates  as  evidenced  by 
a Ca/P  ratio  of  1.583  and  the  report  by  Gee^  that  pyro- 
phosphate was  formed  upon  heating  it  at  600 °C  for  60  hours. 
The  results  reported  by  Rootare  et  a 3^  form  two  distinct 
populations  when  plotted  in  Figure  2.  The  regression 
line  for  the  points  at  25 °C  (the  point  30.20,  19.61 
excluded)  has  a slope  of  -0.778,  which  would  correspond 
to  a ratio  of  Ca/P  = 1.28  in  the  solid  phase.  The  points 
from  treatments  at  40 °C  appear  to  define  two  lines  with 
slopes  -0.309  and  -0.211  which  would  correspond  to  Ca/P 
ratios  of  3024  and  4.74  in  the  solid.  If  equilibrium 
with  respect  to  Ca2 (HP04) (OH) 2 had  been  established,  as 
claimed  by  Rootare  et.  a .]&,  lines  with  slope  of  -0.5  (Ca/P  = 
2)  should  have  been  defined  by  their  data  in  the  potential 
plot.  The  fact  that  their  data  define  two  or  three  Ca/P 
ratios,  none  of  which  correspond  to  their  postulated  surface 
complex,  along  with  the  other  criticisms  listed  above, 
argue  very  strongly  against  the  existence  of  the  surface 
complex  hypothesized  by  them.  LaMef^7  re-examined  the  data 
of  Rootare  et  al°  and  came  to  the  conclusion  that  the  pro- 
duct, (Ca++)  (I1P04  = ) showed  the  least  variation.  However, 
for  the  reasons  given  above,  this  is  not  an  adequate  basis 
for  any  conclusions  regarding  the  nature  of  the  equili- 
brating solid  phase. 


-22- 


In  the  light  of  the  foregoing  results  and  discussion, 
we  conclude  that  our  samples  of  preheated  HA  did  equili- 
brate according  to  the  stoichiometry  of  HA  and  according 
to  the  usual  solubility  product  principles,  that  the 
steam-heated  HA  yielded  a solubility  constant  usable  as 
a frame  of  reference  for  well-crystallized  HA,  and  that 
most  of  the  difficulties  reported  in  the  literature  on 
this  matter  relate  to  contamination  of  the  solid  and  lack 
of  equilibrium  between  the  solution  and  the  solidQ 

7.  APPENDIX 

Ionic  activity  coefficients  were  calculated  from  the 
extension  of  the  Debye-Hiickel  limiting  law. 


log  ft 


AZ3/p, 

1 + BalvA  p, 
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The  values  used^  for  the  temperature-dependent  constants 
A and  B were  0.5092  and  0o3286  x 10s,  respectively;  those 
for  the  parameters  at  were  Ca++,  6 x 10*~8cm;  H2PC>4, 
HPO4,  and  PO4,  4 x 10~8  cm;  H+ , 9 x 10“8 cm;  and  0H“ , 

3.5  x 10“8 cm. 


The  values  used  for  the  three  ionization  constants 
of  E3PO4  were  iq  , 7.108  x 10”3  K2  , 6.31e  x 10“8 

Kg,  4.73  x 10~13  O;  and  for  Kw,  1.012  x 10"1 4 
The  association  constants  used  in  this  paper  are  defined 
by 


[CaHP04 ] 0 

Kx  = . — — _ Ky  = 

(Ca++) (HPOJ) 

with  values17  Kx , 244.0;  and  Ky  , 8o20. 


(CaH2P04+ ) 
(Ca+  + ) (H2 POJ) 
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The  points  in  the  potential  plot  (Figure  2)  corres- 
ponding to  the  data  of  Rootare  et_  al5  at  40 °C  were  cal- 
culated using  the  following  constants:  Kj  , 5099s  x 10"“3  3J> 

Kg,  6 . 607  x 10“8  ^>7  Kg  , 7 . 58s  x 10"1 3 Kw , 20936 

x 10"1 4 ^ ; A (equation  19),  0o5241^>;  and  B,  0.3318 
x 108 
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Fraction  of  Ca  present  as  ion  pairs 
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Solution  Compositions  at  Singular  Points 


Solid  Phases 
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Legend  of  Figures 


Figure  1.  - Solubility  isotherms  of  steam-heated 
HA,  full  circles,  and  air-heated  HA,  open  circles, 
at  25  °C. 

The  two  lighter  curves  are  the  isotherms  of 
CaHP04  •2H2d3v7  , DCPD,  and  CaHP04^,  DCPA.  Compositions 
of  the  calculated  singular  points  A,  B,  C,  and  D are  given 
in  Table  6. 


Figure  2.  - Potential  plot,  calcium  hydroxide  vs. 
phosphoric  acid. 

•,  steam-heated  HA;  ©,  air-heated  HA;  0,  data  from 
Claris ; A. , data  from  Levinskas  and  Neuma^v ; □/  data 

from  Rootare  et.  aj^  at  25  CC;  H,  data  from  Rootare  _et  a 
at  40 °C . 
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